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Abstract 

We review our recent determination of the strong coupling a, from the OPAL data for non-strange hadronic tau 
decays. We find that asirn^) - 0.325 ± 0.018 using fixed-order perturbation theory, and Osirn^) - 0.347 + 0.025 using 
contour-improved perturbation theory. At present, these values supersede any earlier determinations of the strong 
coupling from hadronic tau decays, including those from ALEPH data. 



1. Figure [T] shows a number of recent determina- 
tions of the strong coupling, as(m^), from hadronic tau 
decays. The two values at the top are recent determi- 
nations yj from OPAL data [2J, and have signficantly 
larger errors than all the other determinations shown. 
The reason for these larger errors is twofold: (1) er- 
rors in the non-perturbative part of the sum rules used 
in order to extract CjCmJ) have been systematically un- 
derestimated in earlier works; (2) all other values shown 
in Fig. [T] used ALEPH data, with correlations in which 
the effects of unfolding the spectrum have inadvertently 
been omitted 13] . Here we give a brief overview of our 
analysis f\\^\, comparing it with the standard approach 
used in for instance Ref. fSj. We note that for the values 
shown in Fig. [T] only those of Refs. |[T]|5]|6l are based 
directly on data; all others used estimates for the non- 
perturbative part from Ref. lHJ . 
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2. All tau-based determinations of as{m^) start from 
the sum rule 

r t/s W(s) Pexp(«) = -^-S dzw(z)YlovE(z) 

Jo 27ri Jj-|=,„ 

1 p 

— \ ds w{s) Im nDv(i) , (1) 

in which w{s) is a polynomial weight, Pexp(i) is the in- 
clusive, non-strange spectral function taken from exper- 
iment, nopE(z) is the operator product expansion (OPE) 
expression for the non-strange flavor off-diagonal vac- 
uum polarization with perturbation theory constituting 
the (dimension D -0) dominant part, and 

nDv(z) = nQCD(z) - noPE(z) (2) 

the correction for using nopE(z) instead of the (un- 
known) exact nQCD(z)- It is important to include 
an estimate for the duality-violating (DV) term in 
Eq. ([T]), because duality violations are not small near 
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Figure 1: Recent values for a ,(171^) from hadronic tau decays. For 
references, see Ref. |7|. 

the Minkowski axis in the complex z plane, from which 
the integral on the left-hand side of Eq. ([T]i originates. 
In more physical terms, the OPE does not capture the 
hadronic resonances which are clearly visible in the ex- 
perimental spectrum Pexp('S)- 

While the perturbative part of IIopeCz) is known to 
4th order in asiin^) |8| and the condensate correc- 
tions are expected to give a good estimate for non- 
perturbative corrections for |z| = sq large and away from 
the Minkowski axis, we need to use a model in order to 
estimate the effects of the DV term in Eq. ([T]). Based on 
an extensive study of models for the resonances seen in 
Pexpis), we have used the ansatz (see Ref. 141 and refer- 
ences therein) 

- Im nDv(i) = e-^'-^ sin {a + /3s) . (3) 
n 

3. This very brief overview of the theory (for more 
details and references, see Refs. p^Tl) makes it possible 
to compare the "standard" analysis of Refs. |l2l|5l) with 
ours. In the standard analysis: 

• It is assumed that duality violations can be ne- 
glected if one uses "pinched" weights w(s), which 
always include factors (1 - s/sq)", with n = 2 or 
n = 3. Such weights are thus necessarily polyno- 
mials of at least degree 2. 

• One chooses sq = m^, and 5 pinched weights of 
degrees 3 to 7, generating 5 data points, to which 
one fits 4 parameters, as{inl) and the dimension 4, 
6 and 8 OPE condensates. We note that three of 
the five employed moments have problematic per- 
turbative behavior |9|. 



Figure 2: Vector cliannel fit of Eq. |TJ with w = I using OR\L data. 
CIPT fits are shown in red (dashed) and FOPT in blue (solid). The 
(much flatter) black curves represent the OPE parts of the fits. The 
vertical dashed line indicates the location of . 

in the weight w{s) picks out a term of dimension 
D = 2A:-i-2 in the OPE. This assumption was shown 
to be not self-consistent in Ref. |6|. 

In contrast, in our work: 

• Our main fit uses the data for the vector channel, 
and the weight w - 1, for which no OPE terms 
beyond perturbation theory contribute 

• We let 50 vary over an interval [s„„„, m^], with s„„„ 
determined by the quality and stability of the fit; 
typically, Smin ~ 1-5 GeV^. This provides more 
(correlated) data compared to choosing sq - m^. 

• We include duality violations. Our main fit thus 
has 5 parameters, Q'j(m^) and the 4 parameters of 
Eq. Q for the vector channel. 

• We check consistency using also the axial-channel 
data, and adding other weights up to degree 3 (al- 
ways retaining all required terms in the OPE, i.e., 
to dimension 8.) 

For more discussion of our strategy, we refer to Ref. ID. 

Figure[2]shows the result of a fit to the vector channel 
with weight w - I and = 1.5 GeV^, with all cor- 
relations between the data taken into account. This fit 
leads to (FOPT = fixed-order perturbation theory, CIPT 
- contour-improved perturbation theory) 

as(ml) = 0.307 ±0.018 ±0.004 ± 0.005 (FOPT), (4) 
= 0.322 ± 0.025 ± 0.004 ± 0.005 (CIPT) , 
e-^ = 0.02 ± 0.01 . 



• This assumes that OPE condensates of dimension 
10 through 16 vanish, because a term of order s'' 



'Up to numerically negligible o-, connections to these terms PI. 
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Figure 3: Spectral function predicted by the fit of Fig. |2] compared 
witii OPAL data. CIPT fit.s are shown in red (dashed) and FOPT in 
blue (solid). The (much flatter) black curves represent the OPE parts. 

The value of for this fit is 0.36 per degree of free- 
dom. The first error is the fit error, the second gives an 
indication of the stability with respect to varying .?„„■„, 
and the third estimates the effect of truncating pertur- 
bation theory. From the result for e^^ as well as from 
Figs.[2]and|3]it is clear that duality violations have to be 
taken into account. As mentioned above, and explained 
in detail in Ref. Q, one may attempt to suppress dual- 
ity violations by using pinched weights, but at the price 
of not ignoring terms in the OPE of dimension up to 
16 if the weights of the standard analysis are employed. 
One should then worry, however, whether the OPE con- 
verges to such high order for sq e [.?„„■„, m^]. 

We have carried out a number of consistency checks: 

• Fits with weights w - 1, 1 - (s/sq)^ and (1 - 
s/sofil + 2s /sq) and also including axial data 
lead to results completely consistent with Eq. (j4]|. 
These moments are preferred for their perturbative 
behavior 19]. 

• The first and second Weinberg sum rules, as well 
as the sum rule for the electromagnetic pion mass 
difference are satisfied within errors. 

• Our fits describe the non-strange hadronic tau- 
decay branching fraction Rv+A,ud{-^o) extremely 
well for So above about 1.3 GeV^. 

The values found originally by OPAL are |l2j: 

Usiml) = 0.324 ±0.014 (FOPT), (5) 
= 0.348 + 0.021 (CIPT) . 

We find central values about 0.02 lower on the same 
data, and conclude that OPAL errors were underesti- 
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Figure 4: Ratio of the vector-channel OPAL data using HFAG branch- 
ing fractions and the original OPAL data using 1998 PDG branching 
fractions. The horizontal axis labels the bins. 

mated. In particular, errors due to non-perturbative ef- 
fects are at least as large as the difference between CIPT 
and FOPT 

4. The 1998 OPAL spectral functions were con- 
structed as a sum over all exclusive modes normalized 
with the 1998 PDG values of the branching fractions. In 
Ref 1 1 1, we applied the analysis of Ref. ||4l to rescaled 
OPAL data obtained by instead using current branching 
fractions from HFAG |10|. For the vector channel, the 
rescaling factor is shown in Fig. |4] 

In addition, we also carried out a Markov-chain 
Monte Carlo analysis of the rescaled data in order to 
explore the distribution as a function of the param- 
eters Q'j(m^) and those of Eq. (jsjl in more detail. This 
distribution forms a five-dimensional landscape, since 
it is a function of 5 parameters. Cross sections of the 
full landscape are shown in Fig. [5] in the top figure 
we show the Monte-Carlo generated distribution pro- 
jected onto the ag-x^ plane; in the bottom figure, we 
show the projection onto the ff j-^v plane. 

We see that the;^'- distribution corresponds to a rather 
complicated landscape. The fit appears to allow for 
two different minima, one with 6v ~ 4- and one with 
6v ~ -2. As can be seen in the lower figure, these two 
minima are not well separated, and the value of x^ per 
degree of freedom is reasonable for each minimum. 

While the absolute minimum is at 6v ~ 4, clearly we 
need physical input in order to arrive at a value of a, 
which is more precise than a value in the range 0.27- 
0.35 or so that would follow from the top figure. The 
model developed in Ref. ifTTIl favors 

6v^-\og-^^4.2, (6) 

p 

which leads us to choose the absolute minimum of the 
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Figure 5: Monte-Carlo analysis of the distribution for the fit of 
Eq. ^ with H' = 1 to the vector channel for the rescaled OPAL data. 
V stands for "vector channel." 

distribution as the preferred value, and consider the 
other minimum unphysical. With this choice, we find 

as(ml) = 0.325 + 0.018 (FOPT) , (7) 
= 0.347 + 0.025 (CIPT) . 

While the central values are virtually the same as those 
of Eq. (j5]l, this is purely accidental. Moreover, as before, 
the errors are larger, because of non-perturbative effects. 
Indeed, if we define 6^^ by 

Rv.A^ci = A^c|V,./l' (l + d"''''"- + , (8) 

we find that 

^ -0.004 ±0.012 (FOPT), (9) 
= -0.002 ±0.012 (CIPT) . 

This is to be compared with 6^^ = -0.0059 ± 0.0014 
from the standard analysis ||5][l2l. Our error on 6^^ is 
an order of magnitude larger, because non-perturbative 
eff'ects have been treated systematically in our analysis. 

5. We presented a new analysis of hadronic tau de- 
cays, yielding a new value for the strong coupling, a,,, at 
the tau mass, with a larger error than found in previous 



analyses. Earlier values were all based on the standard 
analysis summarized in Sec. 3, and, moreover, on the in- 
complete ALEPH data. Therefore, we believe that our 
values should be taken as superseding all earlier values 
for as from hadronic tau decays. 

As we saw in Sec. 4, fits to OPAL data are at the very 
edge of what is statistically possible. This is not a flaw 
of the analysis, but appears to be the best one can expect 
based on the OPAL data. In Ref. [1], we investigated 
what would happen with the distribution shown in 
Fig. |5]if the errors are reduced by a factor 2 or 3. We 
found that the degeneracy shown in Fig. |5] disappears. 
Therefore, we expect that a much more stringent test of 
the sum-rule analysis of hadronic tau decays would be 
possible if inclusive spectral functions were made avail- 
able from the BaBar or Belle data. 
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